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ABSTRACT 



Simulations predict that galaxies grow primarily through the accretion of gas that 
has not gone through an accretion shock near the virial radius and that this cold 
gas flows towards the central galaxy along dense filaments and streams. There is, 
however, little observational evidence for the existence of these cold flows. We use 
a large, cosmological, hydrodynamical simulation that has been post-processed with 
radiative transfer to study the contribution of cold flows to the observed z = 3 column 
density distribution of neutral hydrogen, which our simulation reproduces. We find 
that nearly all of the H i absorption arises in gas that has remained colder than 10^'^ K, 
at least while it was extragalactic. In addition, the majority of the H i is falling rapidly 
towards a nearby galaxy, with non-negligible contributions from outflowing and static 
gas. Above a column density of A^hi = 10^^ cm~^, most of the absorbers reside inside 
haloes, but the interstellar medium only dominates for TVhi > 10^^ cm^^. Haloes with 
total mass below 10^'' Mq dominate the absorption for 10^'' < A^hi < 10^^ cm~^, but 
the average halo mass increases sharply for higher column densities. Although very 
little of the Hi in absorbers with A^hi ^ 10^° cm~^ resides inside galaxies, systems 
with A^Hi > 10^^ cm~^ are closely related to star formation: most of their Hi either 
will become part of the interstellar medium before z = 2 or has been ejected from a 
galaxy at z > 3. Cold accretion flows are critical for the success of our simulation in 
reproducing the observed rate of incidence of damped Lyman-a and particularly that 
of Lyman limit systems. We therefore conclude that cold accretion flows exist and 
have already been detected in the form of high column density Hi absorbers. 

Key words: galaxies: evolution - galaxies: formation - galaxies: high-redshift - 
intergalactic medium - cosmology: theory - quasars: absorption lines 



1 INTRODUCTION 

Galaxies grow by accreting gas from their haloes and their 
haloes grow by accreting gas from the intergalactic medium 
(IGM). Feedback from star formation and active galactic nu- 
clei returns a significant fraction of the interstellar medium 
(ISM) to the halo and may even blow it out of the halo into 
the IGM. This cycle of inflow and outflow makes the cir- 
cumgalactic medium and the IGM vital ingredients for our 
understanding of the formation and evolution of galaxies. 

Theoretically we expect that gas accreting onto haloes 
with sufficiently low circular velocities will not shock-heat 
to the virial temperature of the halo, but will instead fiow 
in cold (i.e. T ~ 10* K) and relatively unimpeded. This so- 
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called "cold accretion" will happen if the cooling time of 
virializcd gas is to o short to m a intain a hot, hydrostatic 
halo ( e.g. White fc FrenkI Il99ll : iBirnboim fc D ckcl 20q3; 
iDekel fc Birnboimll2006l '). The existence of such a cold ac- 
cretion mode has been confirmed by simulations, which have 
furthermore demonstrated that cold-mode accretion can also 
be important for haloes sufiiciently massive to contain hot, 
hydrostatic gas. Because gas accretes preferentially along 
the filaments of the cosmic web, the streams of infalling gas 
have relatively high gas densities and correspondingly low 
cooling times. This allows the cold streams to penetrate the 
hot, hydrostatic haloes surrounding massive gala xies, partic- 



ularly at high redshift (e . g. Keres et al. 2005 : Ocvirk et al 
200^; iDekel et al.l l2009l: p^eres et al.l 120091^ iBrooks et al 
12009; 'va n de Voort et all boilallbl: Ivan de Voort fc Schave 



12009; 'va n de Voort et all I2011allbl: Ivan de Voort fc Schs 
I2OII ; Fumagalli et al.ll201ll : lFaucher-Giguere et al.ll201ll ) 
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Although the dilute halo gas that has been shock- 
heated to the virial temperature is routinely detected in 
X-ray observations of clusters and groups of galaxies, and 
has perhaps even b ee n seen around individual gal axies (e.g. 
ICrain et al.] l2010al lbl: lAnderson fc BregmanlbOllT ). there is 
no clear, direct o bservational e viden ce for cold-mode ac- 
cretion. However, iRauch et all (|201ll ) have observed cold, 
filamentary, infalling gas which could be cold-mode accre- 
tion. It has been claimed that the diffuse Lyman-a emis- 
sion detected around some high-redshift g a laxies is due 
to co ld accretion (e.s. IPiikstra fc Loebl[2003 : iGoerdt et al.l 
but both simulations and observations indicate that 
the emission is more likely scattered light from central Hll 
regio n s (e.g. iFurlanetto et al.l I2OO5I: iFaucher-Giguere et al 



2010| : ISteidel et all I2OI0I : iHaves et al.l I2OIII : iRauch et al 
2OIII ). Individual Hi absorbers have also been suggested 



to be evidence for cold accretion base d on their proximity 
to a galaxy and their low metallicity (jRibaudo et al.l [20111 : 
ICiavalisco et al.| l201ll ). but it is difficult to make strong 
statements for an individual gas cloud, particularly since 
Hi near galaxies does not need to be inflowing, even if it 
has a low metallicity. Given that outflows are routinely de- 
tected in the form of blueshifted interstellar ab sorption lines 
in the spectra of star-forming galaxies (e .g. ISteidel et al.l 
I2OIOI : iRubin et al]|2010l : iRakic et al.]|2011al '). one may won- 
der why inflowing gas is not commonly seen in the form of 
redshifted absorption lines. It is, however, quite possible that 
the inflowing material has such small cross-sections that the 
signal is completely swamped by outflowing materi al (e.g. 
iFaucher-Gigufere fc Kereil201ll : lstewart et al.ll2011al ). 

It is challenging to observe cold gas around distant 
galaxies in emission, because it typically has densities that 
are low compared to that of the cold ISM. Absorption 
line measurements are therefore very important, with the 
Lyman-Q line of neutral hydrogen being the most sensi- 
tive probe. Indeed, H I absorption can be detected both 
within and outside of galaxies and at 2 ~ 3 the col- 
umn density distribution has been measured over ten or- 
ders of magnitude in iVH i fe.g.lTvtle JIlQSTlilKim et al.ll2002l: 
Peroux et a"Lll2005l:[OT^ ara ct al. 20 Q7|;|Prochaska fc WoM 
20091 : iNoterdaeme et all I200a : ,Prochaska et all I2OI0I ). By 



correlating the H I absorption in the spectra of background 
quasars with both the transve rse and line of sigh t separa- 
tions from foreground galaxies, iRakic et al.1 (|2011bl ) have re- 
cently presented strong evidence for infall of cold gas on 
scales of ~ 1.4 — 2.0 proper Mpc at « 2.4. 

The low column density (A^'hi < lO'^^ '^ cm~^) material 
is known as the Lyman-a forest and originates mostly in 
the photo-ionized IGM fe.g. iBi et al.l ll992l: 'Cen et al.lll994'; 
iHernauist eral]|l996l : [xheuns et al.1 1998: Schayc 2001a). 
Because these systems are optically thin, thoy arc relatively 
easy to model in simulations. Lyman Limit Systems (LLSs; 



10"'^ < iVHi < 10^"'^ cm"^) and Damped Lyman-a Sys- 
tems (DLAs; Nm > 10^°'^ cm"^) are optically thick to Ly- 
man limit photons. Because the gas comprising these strong 
absorbers is partially shielded from the ambient UV radia- 
tion, it will be more neutral, and will also be colder, than 
if it were optically thin. These high column density systems 
are harder to model, because in order to know which gas 
is self-shielded, one needs to perform a radiative transfer 
calculation. Additionally, at the highest gas densities a frac- 
tion of the hydrogen will be converted into molecules, re- 



ducing the atomic hy drogen fractiorQ (e.g. ISchavd l2001bl : 
iKrumholz et al.ll2009l ). 

The difficulty with interpreting absorption-line studies 
is that they only allow us to study the gas along the line-of- 
sight direction. It is therefore hard to determine what kind of 
objects are causing the absorption. Simulations that repro- 
duce the observed absorption-line statistics can be employed 
to study the three-dimensional distribution of the absorbing 
gas and to guide the interpretation of the observations. 

A lot of progress has been made in understand- 
ing and modelling the H I c o lumn density distribu- 
tion function (e.g. iKatz et al.1 Il996l: ISchavd l2001allb] : 

Ul; 



Zheng fc Miralda^Escudel |2002|: 



McQuinn et al.l 120111 : iFumagalli et al 



Altav et al.1 120111 : 
' 2OIII ). Building 



on previous wo rk, we ask what k ind of gas is causing the 
Hi absorption. lAltav et al.l HoTl) have post-processed the 
reference model of the OWLS suite of cosmological simula- 
tions (see Section [2} with radiative transfer and found that 
the predicted 2 = 3 Hi column density distribution agreed 
with observations over the full range of observed column 
densities. We analyse their simulation, selecting different 
subsets of the gas to see what fraction of the high column 
density absorption is due to what kind of gas. The selections 
are based on the maximum past temperature, membership 
of haloes and the ISM, and the radial velocity towards 
the nearest galaxy. These properties are all important 
to determine whether or not the gas is accreting in the 
cold mode. We find that almost all high column density 
absorption arises in gas that has never been hotter than 
10^'^ K, at least not while it was extragalactic, and that 
resides inside of haloes. Most of this gas is strongly linked to 
star formation because it is either currently in the ISM, has 
been part of the ISM, or will become part of the ISM before 
z = 2. Inflowing gas provides the largest contribution to the 
Hi column density distribution function, but outflowing 
and static gas cannot be neglected. 

In Section [2] we describe the simulation we used. The 
different gas samples are listed and described in Section |3l 
In Section 13] we determine the contribution of the different 
gas samples to the total amount of gas and H I and to the H I 
column density distribution function. We discuss our results 
and conclude in Section [S] 



2 SIMULATIONS 

We use a modifie d version of GADGET-3 (last described 
in ISpringej |2005| ). a smoothed particle hydrodynamics 
(SP H) code that uses the entropy formulation of SPH 
( Spr ingel fc Hernguistl I2OO2I ) , which conserves both energy 
and entropy where appropriate. This work is part of 
the Overwhelming ly Large Simulations (OWLS) project 
l|Schave et al.ll2010l ). which consists of a large number of cos- 
mological simulations with varying (subgrid) physics. Here 
we make use of the so-called 'reference' model, with the dif- 
ference that it was run with updated cosmological parame- 
ters. We have tested runs with different supernova feedback 
and with AGN feedback and we found our conclusions to be 
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ro bust to chang e s in th e model. The model is fully described 
m ISchave et all l|2010l l and we will only summarize its main 
properties here. 

The simulation assumes a ACDM cosmology with pa- 
rameters Qm = I - = 0.272, = 0.0455, h = 0.704, 
as = 0.81, n = 0.9 67. These values are taken from the 
WMAP year 7 data (|Komatsu et al.l 120111 ). In the present 
work, we use the simulation output at redshifts 3 and 2. 

A cubic volume with periodic boundary conditions is 
defined, within which the mass is distributed over 512^ dark 
matter and as many gas particles. The box size (i.e. the 
length of a side of the simulation volume) of the simulation 
used in this work is 25 comoving Mpc. The (initial) par- 
ticle masses for baryons and dark matter are 2.1 x lO'' M0 
and 1.0 X 10^ M0, respectively. The gravitational softening 
length is 2 comoving kpc, i.e. 1/25 of the mean dark mat- 
ter particle separation, with a maximum of 0.5 /i~^proper 
kpc. 

Star formation is mo d elled according to the recipe 
of ISchave fc Dalla Vecchial ()2008l ). A polytropic equation 
of state Ptot oc pgas is imposed for densities exceeding 
nn = 0.1 cm~ , where Ptot is the total pressure and Pgas 
the density of the gas. Gas particles with proper densities 
riH ^ 0.1 cm~'^ and temperatures T 10^ K are moved onto 
this equation of state and can be converted into star parti- 
cles. The star formation rate (SFR) per unit mass depends 
on the gas pr essure and reprod uces the observed Kennicutt- 
Schmidt law (jKennicuttll 1998h by construction. 

Feedback from star formation is impleme nted using the 
prescription of lOalla Vecchia fc Schavell(2008l ). About 40 per 
cent of the energy released by type II supernovae is injected 
locally in kinetic form, while the rest of the energy is as- 
sumed to be lost radiatively. The initial wind velocity is 
600 km s~^ and the initial mass loading is two, meaning 
that, on average, a newly formed star particle kicks twice its 
own mass in neighbouring gas particles. 

The abundances of eleven elements released by mas- 
sive stars and intermediat e mass stars are followed as de- 
scribed in |w lersma et al.l (l2009bh. We assume the stellar 
initial mass function (IMF) of lChabrieil (l2003l). r anging from 
0.1 to 100 Mq. As described in lWiersma et all (|2009al '). ra- 
diative cooling and heating are computed element by ele- 
ment in the prese nce of the cosmic microw ave background 
radiation and the iHaardt fc Madaul (|200ll ) model for the 
UV/X-ray background from galaxies and quasars in the op- 
tically thin limit. 

We apply a self-shielding correction in post-processing 
by using the rad iative transfer calculations described in 
lAltav et al.l l|201ll ) and Altay et al. (in preparation). This 
is an iterative method that uses reverse ray-tracing (i.e. 
shooting rays starting from the gas element) to compute 
the optical depth and to calculate a new neutral fraction 
until the results converge. The method assumes that the ra- 
diation field is dominated by the UV background radiation. 
lAltav et al.l ([2OII) have shown that this provides a good 
match to the observed H I column density distribution func- 
tion. The temperature of the ISM gas in the simulation is 
determined by the imposed effective pressure and therefore 
not realistic, so we set it to 10* K, typical of the warm neu- 
tral medium. Gas with lower densities (tih < 0.1 cm~^) was 
assumed to be optically thin during the simulation, whereas 
this may not be the case because of self-shielding. We use 



the temperatures that were calculated in the presence of the 
UV background to be self-consistent with the simulation, 
even though the actual gas temperatures will be lower in 
the self-shielded case. These choices do not affect the results 
described here. Hydrogen is converted into mole cules with 
a pres cription based on observations bv lBlitz fc Rosolowskvl 
(|2006l) . who provide the molecular fraction as a function of 

the gas pressure. 

Recentlv iKeres et al.l (|201ll ) have shown that the sizes 
of galaxy discs de pend on the com p utatio nal method 
used (but s ee als o iHummels fc BrvarJ (|201ll )). However, 
ISales et al.l l|2010l ) have shown that changes in the feed- 
back model lead to dramatic variations in the morphology 
of discs, as well as to variati ons of up to an ord e r of m agni- 
tude in disc masses. Indeed. [Scannapieco et ahl (|201ll l have 
shown explicitly that variations in the feedback prescription 
are much more important than the numerical scheme used 
to solve the hydrodynamics. To assess the robustness of our 
conclustions, we have redone our analysis using simulations 
with different feedback prescriptions, but do not show the 
results. The simulations reproduce the observed Hi column 
density distribution function independent of the feedback 
scheme (Altay et al. in preparation). Although the relative 
contributions of the different gas selections to the Hi ab- 
sorption vary somewhat between different simulations, the 
general conclusions are unchanged. We therefore conclude 
that the uncertainties associated with the feedback and the 
numerical method are unlikely to be important for our con- 
clusions. 



3 GAS SAMPLES 

We compute the z = 3 contribution of gas with different 
properties to the H I column density distribution function. 
The gas selections that we consider are listed in Table [1] 
They are based on the maximum past temperature, Tmax, 
membership of haloes, the peculiar velocity towards the 
nearest halo, Wrad, and on whether the gas is part of the ISM. 
These properties are all important to determine whether or 
not this gas is accreting in the cold mode. Cold-mode gas 
should have a low temperature at least until it becomes part 
of the ISM. Gas has only been accreted if it is inside a halo. 
It is likely to accrete onto a halo or its central galaxy if its 
radial peculiar velocity towards the halo centre is high, but 
note that gas orbiting around the central galaxy can also be 
built up by col d flows and could be viewed as evidence for 
cold accretion (jStewart et al.l [ioilbh . The gas is certainly 
accreting onto a galaxy if it enters the ISM of a galaxy and 
becomes star forming in the near future. Below, we describe 
the selections in detail. 

The maximum past temperature, Tmax, is the maximum 
value that the temperature of a gas element has reached 
over the entire past simulation history. High time resolution 
is achieved by storing a separate variable during the simu- 
lation which is updated every time-step if the gas reaches 
a higher temperature than the previous value of Tmax. The 
artificial temperature that a gas particle has while it is on 
the imposed equation of state, and thus star forming, is ig- 
nored. Because we are interested in the temperature history 
before the gas accretes onto a galaxy, this is consistent with 
our aims. Cold-mode gas does not experience a shock at 
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Table 1. List of selection criteria for the various gas samples in the same order as in Figures [T] and |4l 



# 


name 


selection 


1. 


« 1 ] 

all 


all gas 


2. 


Tmax<5p5 


gas with Tmax ^ 10^-^ K 


3. 


inhalo 


nalo gas (witn iVi^alo ^ J^U MqJ 


4. 


inmainhalo 


gas m mam haloes (with M^^j^ jiaio ^ 10 Mq) 


5. 


insubhalo 


gas in satellites (with Mgatellite ^ 10^ 1^©) 


6. 


intiaLoytolu 


gas in haloes with lU ^ J^halo < lU NLq 


7. 


intiaLol Utol 1 


gas m haloes with iU ^ ^"halo 


8. 


mnalolltolz 


gas in haloes with lU ^ ^^ihalo < lU Mq 


y. 


tntialo 1 zto 1 J 


gas in haloes with lu -^^halo -l'-' ■'■^^^0 


10. 


inflowing 


gas that is moving towards the nearest (nearest in units of /2vir) halo centre with a 






peculiar velocity greater than a quarter of its circular velocity 


11. 


static 


gas that is neither mnowmg (#10) nor outnowing (#12) 


12. 


outflowing 


j.lJ_" ■ £ J.1 J./ J." "j-TONll j_ "j-l 

gas that IS moving away from the nearest (nearest m units ol -Kvir) halo centre with 






a peculiar velocity greater than a quarter of its circular velocity 


13. 


inhalo^inflowing 


halo gas that is moving towards the halo centre with a peculiar velocity greater than 






a quarter of its circular velocity 


14. 


inhalo^static 


halo gas that is neither inflowing (#10) nor outflowing (#12) 


15. 


inhalo^ outflowing 


halo gas that is moving away from the halo centre with a peculiar velocity greater 






than a quarter of its circular velocity 


16. 


inISM 


gas which is in the ISM at 2: = 3 


17. 


pastISM 


gas which is not in the ISM at 2; = 3, but was at some redshift z > 3 


18. 


futureISM 


gas which is not in the ISM at 2 = 3, but will be before 2 = 2. 


19. 


futurefirstISM 


gas which is not in the ISM at 2 ^ 3, but will be before z = 2. 



the virial radius and reaches a maximum temperature of 
~ 10^ K just before accretion onto a halo due to heating by 
the UV background ijvan de Voort et"alll2011al ). We there- 
fore select cold-mode gas by using Tmax ^ 10^'® K (sample 
Tmax<5p5 of Table While hot- mode gas has by defini- 
tion been hotter th an 10^'^ K, its current tempe rature can 
be much lower (e.g. Ivan de Voort fc Schavell201ll ). 

We note that Tmax may underestimate the true max- 
imum past temperature. In SPH simulations a shock is 
smeared out o ver a few smoothing length s, leading to in- 
shock cooling (|Hutchings fc Thomad [2OO0I ). If the cooling 
time is of the order of or smaller than the time step, 
then the maj ci mum temperature will be underestimated. 
ICreasev et alj (|201ll l have shown that a particle mass of 
10^ Mq is sufficient to avoid numerical overcooling of accre- 
tion shocks onto haloes, which is very close to the resolution 
used in this simulation. Interpolating their results indicates 
that our mass resolution is indeed sufficient. Even with infi- 
nite resolution, the post-shock temperatures may, however, 
not be well defined. Because they have different masses, elec- 
trons and protons will have different temperatures in the 
post-shock gas and it may take some time before they equi- 
librate through collisions or plasma effects. We have ignored 
this complication. Another effect, which was also not in- 
cluded in our simulation, is that shocks may be preceded 
by the radiation from the shock, which may affect the tem- 
peratur e evolution. Ignoring these issues. Ivan de Voort et aP 
l|2011al ) showed that the distribution of Tmax is bimodal and 
that a cut at Tmax ~ 10^'^ K naturally divides the gas into 
filamentary cold- and diffuse hot-mode accretion and that it 



produces similar results as studies based on adaptive mesh 
refinement simulations. 

Haloes are found by using a Friends-of-Friends (FoF) 
algorithm. If the separation between two dark matter par- 
ticles is less than 20 per cent of the average separation (the 
linking length b — 0.2), they are placed in the same group. 
Assuming a radial density profile p(r) oc r~^, corresponding 
to a flat rotation curve, such a group has an averag e overden- 
sity of (phaio>/(p> ~ 180 (e.g. lLacev fc ColdllQM ). which is 
close to the value for a virialized object predicted by the top- 
hat spherical collapse model. Baryonic particles are placed 
in a halo if their nearest dark matter neighbour is part of 
the halo. Haloes can contain both bound and unbound par- 
ticles. The selections are made by using haloes with total 
mass Mhaio ^ 10® M0 , corresponding to approximately 100 
dark matter particles. A gas particle is considered to be in a 
halo if it is a member of a Friends-of-Friends group {inhalo). 

W e use SUBFIND ISpringel et all I2OOII : iDolag et al.l 
I2OO9I ) on the FoF output to identify gravitationally bound 
subhaloes. The main halo is the most massive subhalo in 
a FoF halo and the remaining subhaloes are classified as 
satellites. Because only bound particles are included in this 
definition, some unbound particles are not attached to any 
subhalo, although they are part of the FoF halo. The gas 
is in a main halo if it is a member of the most massive 
subhalo, with A/main halo ^ 10^ {iumainhalo) and the 
gas is in a satellite if it is a member of a subhalo, with 
Afsateiiite ^ 10^ M© , that is uot a main halo (insubhalo). 

The halo gas (inhalo) is subdivided into four mass bins, 
lO'-* ^ A/halo < 10^" Mq (mhalo9tol0), 10^° «: Afhaio < 
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10" M© (inhalolOtoll), lO" ^ Mhaio < 10^^ M© {m- 
halolltol2), and 10^^ Mhaio < 10^^ Mq {mhalol2tol3), 
containing 38 504, 3626, 275, and 17 haloes, respectively. 

The radial peculiar velocity of every gas particle is cal- 
culated with respect to the nearest halo. The position (and 
the centre) of the halo is defined as the position of its most 
bound particle. The virial radius, i?vir, is the radius around 
this centre which encloses an average density that agrees 
with the prediction from top- hat spherical collapse ca lcula- 
tions in a ACDM cosmology l|Brvan fc Normanlll998l '). Dis- 
tances towards a halo, d, are normalized by Rvii of that halo. 
The nearest halo is the one for which d/R^ir is minimum. 
The halo's peculiar velocity is defined to be the centre of 
mass velocity of its main subhalo and is subtracted from 
the gas peculiar velocity to determine the relative radial 
peculiar velocity. If the gas is located inside a halo, its ra- 
dial velocity is calculated with respect to the centre of that 
halo, even if it is also in a satellite. The selections are made 
by comparing the radial velocity of the gas to a quarter of 
the circular velocity, Wcirc = \/ GMvir/-Rvir, of the nearest 
halo, where Mvir is the total mass inside -Rvir. The gas is 
inflowing when Urad ^ —0.25110110 (inflowing), static when 
— 0.25ucirc < Viad < 0.25i'circ (static) , and outflowing when 
irad 0.25«circ (outflowing) . Notc that gas near or inside a 
small halo embedded in a fllament may be classified as out- 
flowing with respect to its nearest halo, whereas it might be 
inflowing with respect to a much larger halo that is being 
fed by this filament. It is hard to unambiguously identify 
this gas, so we just caution the reader that the amount of 
inflowing gas may be an underestimate. This would only 
strengthen our conclusions. 

The definition of the nearest halo is somewhat arbitrary 
and the gas may be many virial radii away and unrelated 
to the nearest halo. This is why we make the same radial 
velocity cuts for all gas (samples 10-12 of Table [T]) also for 
gas that is already inside a (FoF) halo, so that we know 
it has already been accreted by the halo (samples 13-15 of 
Table [T]). Halo gas which is inflowing is in the process of 
accreting onto the central galaxy [inhaloAnflowing). 

In our simulation gas becomes star forming when its 
density exceeds nn ~ 0.1 cm~^ while its temperature T ^ 
10^ K. Because of supernov a feedback, the gas can leave 
the ISM and re-enter it fe.g. lOppenheimer et al.ll2010l \ We 
will call this recycled gas even though it may not have been 
part of a star before it was ejected from the ISM. We keep 
track of the time a gas particle last reached or left the ISM. 
We select all the gas that is star forming at z = 3, i.e. the 
ISM (inlSM), star forming before but not at z = 3, i.e. 
recycled gas (pastlSM), star forming before z — 2 but not 
at z = 3, i.e. gas that will accrete onto a galaxy before z — 2 
(futurelSM), or star forming before z = 2 but not at 2 ^ 3, 
i.e. gas that will accrete onto the galaxy before z = 2 for 
the first time (futurefirstlSM). The time elapsed between 
z — 3 and z — 2 is about 1.2 Gyr, which is similar to the gas 
consumption time scale implied by the observed Kennicutt- 
Schmidt star formation law for typical ISM densities and 
indicates whether the gas will accrete onto a galaxy in the 
near future. 



4 RESULTS 

4.1 Gas and Hi fractions 

Many of the gas samples listed in Table [T] contain a much 
larger fraction of the total Hi mass than of the total gas 
mass. This is illustrated in Figure [Tl which shows the frac- 
tion of the total amount of gas (left panel) and the fraction 
of the total amount of neutral hydrogen (right panel) ac- 
counted for by the samples listed in Table [T] The density 
parameter of all the gas Slgas ~ 0.15787 and of neutral hy- 
drogen r^Hi = 0.00178, so at z = 3 only about one per cent 
of the gas mass consists of neural hydrogen. Selections based 
on halo membership (samples 3-9 and 13-15) and on partic- 
ipation in star formation (samples 16-19) account for little 
of the gas mass, but a significant part of the Hi mass. 

The maximum past temperature selection ( Tmax<5p5) 
contains almost all gas in the simulation. This is not so sur- 
prising, because most of the mass in the universe is in the 
IGM and the UV background can only heat this gas to tem- 
peratures < 10^ K. The same selection also contains almost 
all neutral hydrogen gas. Gas that shocked to T > 10^'^ K 
while it was extragalactic, either through accretion shocks 
or shocks caused by galactic winds, accounts for only about 
10 per cent of both the gas and Hi. This already tells us 
that most gas and most neutral hydrogen is cold-mode gas. 

Note that this is not a trivial result. While it is true 
that gas with temperatures above 10^'^ K will be collision- 
ally ionized and therefore difficult to see in Hi absorption, 
we distinguish cold- and hot-mode gas based on the maxi- 
mum past temperature rather than the current temperature. 
Gas parcels accreted in the hot mode can cool down, reach 
temperatures similar to those of cold-mode gas, and show 
up in Hi absorption. 

Haloes contain less than 10 per cent of the gas, but 
more than 90 per cent of the neutral hydrogen (inhalo). Gas 
outside haloes contributes little to the Hi content, because 
the neutral fraction increases with density and the highest 
densities are found inside haloes. Combining this with the 
finding that sample Tmax<5p5 accounts for nearly 90 per 
cent of the neutral hydrogen, proves that most H I absorp- 
tion arises from gas that has been accreted onto haloes in 
the cold-mode. All main haloes contain seven times more 
gas and four times more Hi than all satellites (compare in- 
mainhalo and insubhalo). 

Haloes from the entire mass range contain a signif- 
icant part of the Hi mass (samples 6-9). Haloes in the 
highest mass bin are very rare and their total contribution 
is the smallest (15 per cent; inhalo 12to 13). Haloes below 
Afhaio = 10^ M0, and therefore below our resolution limit, 
are not expected to contain a large amount of H I, because 
we already account for the vast majority of H I with higher 
mass haloes [inhalo) a nd our simulations reproduce the ob- 
served Hi distribution (|Altav et al.ll201ll ). 

Approximately 40 per cent of the gas and half of the 
neutral hydrogen is inflowing faster than a quarter of the 
circular velocity of the nearest (nearest in units of i?vir) halo 
[inflowing). The contribution is even larger for all inflowing 
gas, without a minimum radial velocity threshold (about 65 
per cent; not shown). Outflows always account for less than 
25 per cent of the gas and Hi mass [outflowing). The con- 
tribution of outflowing gas and of static gas (i.e. gas moving 
either in or out at velocities below 0.25i;circ; static) is lower 



© 2011 RAS, MNRAS 000.[TlfT2l 




Figure 1. Fraction of the total amount of gas (left) and Hi (right) at 2 = 3 accounted for by the selections listed in Table [T] Many 
samples (#3—9 and #13-19) contain a much larger fraction of the Hi mass than of the total gas mass. 



for H I than for all gas, whereas the contribution of inflowing 
gas is higher. 

Within haloes, the importance of inflows {in- 
halo -inflowing), with respect to outflows (inhalo-outflowing) 
and static gas (inhalostatic) , increases slightly. This is not 
visible in Figure [Tl because the bars are too small. The H I 
inside haloes is distributed over the inflowing, static, and 
outflowing components in the same way as all H I (i.e. both 
inside and outside of haloes), because haloes account for 
nearly all the Hi. Thus, about half of the Hi can be ac- 
counted for by gas inside the virial radius that is falling 
towards the halo centre with a velocity greater than 25 per 
cent of the halo's circular velocity. 

Very little gas is contained in the ISM of galaxies, but 
almost half of the Hi mass is (inlSM). The amount of gas 
that has been ejected from galaxies at 2; > 3 (pastlSM) is 
greater than the amount that is in the ISM at 2: = 3 (inlSM), 
but the situation is reversed for Hi. Similarly, more gas has 
been ejected from galaxies at 2 > 3 (pastlSM) than will 
accrete onto galaxies at 2 ^ z < 3 (futurelSM), but it is the 
other way around for neutral hydrogen. Taken together, gas 
that is currently in the ISM and gas that will enter the ISM 
within 1.2 Gyr {inISM + futureISM) account for most of the 
Hi in the universe. Gas/Hl that will accrete onto a galaxy 
at 2 2 < 3 for the first time (futurefirstlSM) accounts for 
most/less than half of the total amount of gas/Hl that will 
accrete at 2 ^ 2 < 3 (futureISM). 

4.2 Spatial distribution: A visual impression 

To gain insight into the spatial distribution of the neutral 
hydrogen. Figure [2] shows the Hi column densities at 2 = 3 
in a 2 by 2 comoving Mpc region centred on a lO^'^'^ M0 
halo, one of the most massive objects in the simulation. The 
Figure is meant to illustrate the spatial distribution of the 
H I corresponding to some of the samples listed in Table [T] 
Some gas samples trace only a subset of the high column 



density sightlines, while others trace the full morphology of 
the Hi. 

In the top row, the column densities are shown when in- 
cluding all gas (left), only gas with Tmax ^ 10^'^ K (middle), 
and only gas inside haloes (right). The vast majority of neu- 
tral gas inside and outside haloes satisfies Tmax 5; 10^'^ K, 
so even for the 10^^ * Mq halo with Tvir = lO'*'^ K, most 
of the Hi is contributed by gas that has Tmax ^ 10^'^ and 
that has thus not gone through an accretion shock near the 
virial radius. When we exclude gas outside haloes (top right 
panel), we lose the filamentary material connecting different 
haloes, as well as most of the A''hi < 10^^ cm~^ sightlines. 

The middle row shows the column density when includ- 
ing only gas that is inflowing towards the nearest (nearest 
in units of Rvii) halo faster than a quarter of the circular 
velocity of that halo (left), only gas outflowing faster than 
a quarter of the circular velocity (right), and the remaining 
and hence "static" gas (middle). The morphologies of the H I 
gas are very different. The inflowing gas traces most of the 
extended and filamentary structure of the total H I distribu- 
tion, while the static and outfiowing gas are concentrated in 
the centres of the haloes. 

The images shown in the bottom row include only gas 
that is not part of the ISM at 2 = 3. Additionally, the panels 
include only gas that was part of the ISM for some redshift 
2 > 3 (i.e. ejected gas; left panel), only gas that becomes part 
of the ISM at 2 5C 2 < 3 (middle), and only gas that joins 
the ISM at 2 ^ 2 < 3 for the first time (right). The ejected 
gas that is not currently part of the ISM (bottom left panel) 
is located relatively close to the centres of haloes, where star 
formation is taking place, but it is more spread out than the 
actual ISM (not shown). Comparing it to the middle row, 
with inflowing (left), static (centre), and outflowing (right) 
gas, we immediately see that some of the ejected gas is al- 
ready inflowing and thus recycling. Note that it may have 
been ejected from a different halo than the one it is in now. 
Most of the lower column density filamentary gas, especially 
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Figure 2. Hi column density images for a cubic region of 2 comoving h ^ Mpc on a side centred on a 1012-4 M0 lialo at ^ = 3, one 
of the most massive haloes in this simulation. For reference, the virial radius of this halo is 268 comoving kpc and indicated with 
black circles. For all but the first panel, the Hi column densities were computed including only gas selected to be part of the samples 
indicated in each panel and listed in Table [T] Some gas samples trace the full H I morphology, whereas others only trace a subset of the 
high column density sightlines. 



the gas outside the halo, is missing. Even though ejected gas 
contributes significantly to A'hi ^ 10^* cm"^ gas, it is unim- 
portant for lower column densities. 

If we include all the gas that will reach the ISM between 
2 ^ z < 3, i.e. within 1.2 Gyr from the redshift at which we 
are studying the H I column densities (bottom middle panel), 
we leave the morphology of the A^'hi > 10^'' cm"^ sightlines 
intact, including the filaments in between haloes, although 
they become somewhat weaker. Most, but not all, Hi gas in 
filaments close to a halo will accrete onto a galaxy before 
z = 2. 

The last panel shows a subset of the gas in the previous 



panel, because it only includes gas that reaches the ISM be- 
tween 2 ^ z < 3 for the first time. The general morphology 
of H I is the same, apart from some extra holes in the cen- 
tral region of the halo, which are dominated by ejected gas. 
The column densities above 10^* cm"^ are lower for a given 
pixel. 

We conclude that the large-scale H I filaments that feed 
massive haloes consist mostly of gas that has never gone 
through a virial shock, that is inflowing and that will join 
the ISM, and hence participate in star formation, in the near 
future. This gas typically has H I column densities similar to 
those of LLSs or higher for embedded clumps. Within haloes. 
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Figure 3. Predicted Hi column density distribution function at 
2 = 3 (black curve) and a compilation of observation s around 
the same redshift (data points : Pcroux ct al. 2005; O'Me ara et alj 
I2OO7I: iNoterdaeme et al Il2009l: IProch aska ct al. 2010) andpower- 
law constraints (dashed curve: IPro chaska & Wolfe 2009). The 
black curve is obtained by collapsing the full three-dimensional 
H I distri bution onto a pla ne. The observations are matched well 
(see also lAltav et al.ll2011^ . 



the infalling cold-mode gas has higher Hi column densities 
and tends to have a more clumpy morphology than the H I 
in filaments outside haloes. 



4.3 The column density distribution function 

The Hi column density distribution function is defined as 
the number of absorption lines Af, per unit column density 
dA'ni, per unit absorption distance dX. 



dz 



(1) 



dNuidX dNn.dzdX 

Absorbers with a fixed proper size and a constant comov- 
ing number density are distributed uniformly, in a statisti- 
cal sense, per unit absor ption distance along a line-of-sight 
(jSahcall &: Peeble^l 19691 ). The absorption distance is related 
to the redshift path dz as dX/dz = Ho{l + z)'^ /H{z), where 
H{z) is the Hubble parameter. 

The predicted Hi column density distribution func- 
tion is shown in Figure [3] as the black curve. It was 
calculated by projecting the full three-dimensional gas 
distribution onto a two-dimensional grid with 16 384^ 
pixels, which is sufficient to achieve convergence. Ob- 
servationall y determined p ower-law slopes at A^hi < 

and data points at 



10 



A^Hi > 10"" cm 



Prochaska fc W olfe 2009,^ 

^ '' Perouxet al.ll20ol: fO'Meara et al.ll2007l : 



INoterdaeme et al] l2009l : Prochaska et al. 2010h have been 
converted to the cosmology assumed in our simulation and 
are shown as red, dashed lines and red data points, respec- 
tively. The observations of the Hi column density distri- 
bution function are m atched well, also a t column densities 
lower than shown here (|Altav et al.ll201ll ). The match would 
have been even better if we had correct ed the gas temper - 
atures for the effect of self-shielding (see lAltav et al.ll201ll ). 
Including the radiat ion from lo cal sources m ight worsen the 
agreement, alt hough iNagamine et al.l (|2010l ) have shown this 
effect to be unimportant. 



For TVhi > 10"^ cm ^ nearly all sightlines are domi- 
nated bj;_a^_smgle_ateor£fensj;^tOT are very 
rare jProchaska et al.ll20ld : lAltav et alll201ll ). For column 
densities < 10^® cm~'^ our method of projecting the entire 
25 comoving Mpc simulation box does not recover the 
true column density distribution due to t he projection of un - 
related lower column density absorbers. lAltav et al.l ()201ll ) 
therefore determine the column density distribution in this 
regime by decomposing Hi absorption spectra into Voigt 
profiles. Above A^hi ~ 10^* cm~^ the column density dis- 
tribution f unction flattens, because self-shielding becomes 
important jKatz et al.| [ l996l : IZheng fc Miralda-Escudil2002l : 



lAltav et al]l201lh . It steepens again at A'hi ^ 10 

becau se the neutral hydrogen fraction saturates (jAltav et al.l 
lloij), and at Ak^_> 10^ m~^ due to the forma- 



tion of molecules |Schayg 2001bl: Zwaan fc Prochaskall2006l : 
iKrumholz et al" I I2OO9I : INoterdaeme et ai]|2009D . 

We can quantify the contribution of different gas sam- 
ples by calculating the mean fraction of the H I column den- 
sity that is due to each gas selection as a function of the 
total H I column density. This is shown in Figure (Jj for all 
the samples listed in Table [1] Note that the mean fraction 
would be 0.5 for a column density A'^^, if the selection in- 
cludes all A^hj absorbers, but only half of their Hi. However, 
a mean fraction of 0.5 would also be obtained if the selection 
includes half of the absorbers with column density A^hu but 
all the H i of the selected absorbers. Samples based on cuts 
in halo mass are best described by the second case, but most 
other samples combine the two types of possibilities. 

Figure |4] shows that most of the H I absorption is due 
to cold-mode gas that is ffowing towards the centre of the 
nearest halo. Above A'hi = 10^'^ cm~'^ this gas is primarily 
inside haloes. Most of the H I gas at these column densities is 
currently in the ISM or will become part of the ISM within 
the next 1.2 Gyr. We discuss this in more detail below. 

For absorbers with A^hi < 10^" cm~^ more than 90 
per cent of the Hi column density is due to gas that has 
never gone through a virial shock, i.e. cold-mode gas (top left 
panel). For higher column densities, which arise mostly in 
the ISM of galax;ies residing in haloes with Mhaio > 10^^ M0, 
this cold-mode fraction drops to 70 per cent. We do not know 
what caused the remaining gas to be heated above 10^'^ K 
(note that the same gas had a much smaller neutral fraction 
at the time when it was this hot). It could have been a virial 
shock, an accretion shock at a smaller radius, or a shock 
associated with a galactic wind. Note that the gas that has 
Tmax > 10^'^ K because it has been shocked in an outffow, 
may have been accreted in the cold mode. 

Most high column density gas resides in haloes, but the 
halo fraction drops below 50 per cent for A'hi < lO'^^ cm~^ 
(top middle panel). Combining this with the fact that most 
gas has r,„ax < lO'-^ K, we conclude that most LLSs and 
most DLAs arise from gas inside haloes that has been ac- 
creted in the cold mode. Most of the Hi absorption occurs 
in main haloes. Satellites contribute less than 20 per cent for 
all column densities, but their contribution increases slightly 
with the H I column density. 

Very low- mass haloes (10^ < Mhaio < lO'" Mq) are 
most abundant and therefore contain the largest amount 
of Hi below A'hi = lO^'^ cm~^, accounting for about 40 
per cent of the total absorption (top right panel). For 
A'hi > lO'^" cm"'^, a single very low- mass halo is simply 
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not large enough to provide a significant cross section. The 
contribution from haloes with Mhaio > 10^^ M© is dom- 
inant and steeply increasing at A'hi > lO'^^ cm~^. The 
fact that low- (high-) mass haloes dominate low (high) col- 
umn d ensity DLAs agrees qualitatively with'Pontz en et alj 
l|2008l ). although they found that haloes more massive than 
10" '^ Mq only dominate for A^hi > lO^^'^ cm'^ Unre- 
solved haloes (Mhaio < 10® M©) cannot dominate the ab- 
sorption for A^Hi > 10^^ cm~^, because most of it is already 
accounted for by resolved haloes. 

We also split the gas into several SFR bins, but do not 
show the result. We found that for Ahi < 10^^ cm~^, the 
largest contribution (about 30 per cent), comes from haloes 
with 0.01 < SFR < 0.1 M0/yr. For higher column densities 
objects with SFR > 10 Mo/yr are most important, account- 
ing for 30-65 per cent of the Hi absorption, depending on 
the Hi column density. 

At all column densities, most Hi gas is inflowing, but 
there are important and roughly equal contributions from 
static and outflowing gas (bottom left panel). The contribu- 
tions from inflowing, static, and outflowing gas are nearly 
constant (about 55, 30, and 15 per cent, respectively) up to 
Ah I = 10^^ cm~^, at which point the contribution of static 
gas increases relative to those of both in- and outflows. If 
we separate the gas into inflowing or outflowing without any 
minimum radial velocity threshold (not shown), then the 
ratio of inflowing to outflowing remains approximately the 
same, with 65 per cent (35 per cent) of the Hi absorption 
due to inflowing (outflowing) gas. 

Above A'hi = 10^* cm~^, the contribution to Hi ab- 
sorption from gas that is in haloes and inflowing, static, or 
outflowing (bottom middle panel) is almost equally large as 
that from all gas that is inflowing, static, or outflowing (bot- 
tom left panel). Below Ahi ~ 10^* cm~^, their contribution 
is lower, because many absorbers reside outside of haloes. 

The fraction of H I absorption due to gas inside galaxies, 
i.e. the ISM, is only significant for DLAs and dominates at 
A^Hi 10^^ cm~^ (bottom right panel). Two-thirds of the 
10"-5 < Ah, < 10^^ cm material is not participating in 
star formation yet, but will do so within the next 1.2 Gyr. 
About half of this gas will accrete onto a galaxy for the 
first time. Gas previously ejected from galaxies contributes 
a comparable amount of absorption as gas accreting onto 
galaxies for the first time. 



5 DISCUSSION AND CONCLUSIONS 

After post-processing with radiative transfer, the OWLS ref- 
erence simulation matches the observed z = 3 Hi column 
density distributio n over ten orders of magnitude in Ahi 
l|Altav et al.ll201ll ). This success gives us confidence that we 
can use this simulation to study the relation between cold 
accretion and high column density H I absorbers. 

Like other simulations, our simulation also shows 
cold streams and filamen ts inside and outside of haloes 
llvan de Voort et"ai]|2011al '). Gas in many of these streams 
does not go through an accretion shock near the 
virial radius and, because it is able to accrete onto 
a galaxy much more efficiently than the hot, dif- 
fuse gas, i t is v ital for fuell i ng st a ,r formation (e.g . 
iKeres et al.l l2005l : lOcvirk et al] l2008l : iKeres et all l2009l : 



van de Voort et al.l l2011allb[ : Ivan de Voort fc Schavj I2OIII : 
Faucher-Giguere et al.ll201ll ). Here, we demonstrated that 



our simulation would not have matched the observed H I col- 
umn density distribution without cold accretion. This fact 
alone provides evidence that cold accretion has already been 
observed. 

Motivated by the theoretical framework of cold, fila- 
mentary accretion, we have investigated the importance of 
gas that satisfies certain selection criteria for column den- 
sities above A^hi ~ 10^^ cm~^ at 2; = 3. It is difficult to 
check the morphology of all the H I in the simulation, but 
from images (see also Figure [2]) we do know that a lot of the 
accreting cold-mode gas is filamentary. However, the highest 
column density gas tends to be clumpy. Even though much 
of this gas is infalling and has never gone through a virial 
shock, it would not all be classified as streams based on its 
morphology. 

We will first discuss the results separately for three col- 
umn density regimes: i) the Lyman-a forest; ii) LLSs and low 
column density DLAs, and iii) high column density DLAs. 
We will then discuss our results in the context of the cold- 
mode accretion paradigm, compare the m to recent theoret- 
ical results from lFumagalli et al.l |20l3), and conclude. 

(i) Nearly all Hi absorbing gas in the denser part of 
the Lyman-Q forest (Ahi ~ 10^®~^^ cm~^) has never gone 
through a virial shock (Tnmx < 10^'^ K) and most of it is 
outside of haloes (the same is true for the low column den- 
sity forest, but we did not show this here). Even though 60 
per cent is fiowing towards the nearest halo (nearest in units 
of -Rvir) with velocities greater than a quarter of its circular 
velocity, only 15-40 per cent will participate in star forma- 
tion before z = 2, i.e. within 1.2 Gyr. Most of the Lyman-Q 
forest gas that will accrete onto a galaxy before z = 2 will 
do so for the first time. 

(ii) Nearly all Hi absorbing gas in LLSs (10''' < Ah, < 
10^°-^ cm"^) and in low column density DLAs (10^°'^ < 
Ah I < lO'^' cm~^) has never gone through a virial shock. 
About 60-95 per cent of the Hi in LLSs is inside haloes 
with larger contributions from lower-mass haloes: 30-40 per 
cent from 10®"" Mq, 20 per cent from 10^"" Mq, 10-25 



per cent from 10 



M0, and less than 15 per cent from 



■|^gi2--i3 ]y|-^ halocs. Evcu iusidc haloes, almost half of the 
neutral gas is infiowing with velocities greater than 25 per 
cent of the circular velocity. For LLSs, more than 90 per 
cent of the H I is currently outside of galaxies, but 40-70 per 
cent will accrete onto a galaxy within 1.2 Gyr and just over 
30 per cent will do so for the first time. The gas reservoir 
that is feeding galaxies is thus, at least in part, associated 
with LLSs. A significant fraction of the gas in LLSs and low 
column density DLAs has previously already been inside a 
galaxy. This agrees with the finding that the re-accretion of 
gas that has been ejected is impo rtant for the build-up of 
galaxies (|Oppenheimer et al.ll2010l ). 

(iii) About 10-30 per cent of Hi absorbing gas in high 
column density DLAs (Ahi > 10^' cm~'^) was heated above 
■|^g5.5 ^j^g past. This gas has either been accreted in the 
hot mode or it has been shock-heated by a galactic wind. 
Depending on the column density, 55-95 per cent of the 
high- Ah I DLA gas is contained in the ISM of galaxies and 
50-70 per cent of the H I is inside haloes more massive than 
10^' M0. Outfiowing gas contributes about 15 per cent of 
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Figure 4. Mean fraction of the H I column density at 2 = 3 contributed by the gas which fulfils the selection criteria listed in Table [T] 
as a function of the total Hi column density. Almost all of the contributing gas has Tmax < 10^'^ K although hot-mode gas becomes 
significant, accounting for up to 30 per cent of the absorption, above lO^'^ cm~^. More than 80 per cent of A^hi > 10^* cm~-^ gas is 
located inside (main) haloes. About 55 per cent of the H I is inflowing, but the contributions from static and outflowing gas are about 
30 and 15 per cent, respectively, and therefore not negligible. The contribution from gas accreting onto galaxies, i.e. gas that will join 
the ISM in the near future, is dominant for lO^'^ — 10^^ cm~^, accounting for up to 70 per cent. Strong DLAs (Afpii > 10^^ cm~^) 
are dominated by ISM gas (contributing more than 60 per cent to the absorption) and the Lyman-a forest by gas that will remain 
intergalactic down to at least z = 2 (contributing more than 60 per cent to the absorption). 



the Hi. Inflowing gas provides the largest contribution to 
the Hi of high column density DLAs, but the static gas 
component approaches the importance of the inflowing gas 
at the highest column densities. 

The simplest way of defining cold-mode accretion is 
purely based on its maximum past temperature. Gas with 
temperatures below 10^'^ K is able to cool efficiently, be- 
cause the cool i ng fun ction peaks at T « ]^q5-5.5 (e.g. 
|w iersma et al.l l2009al ). Above this temperature, cooling 
times become longer. We define cold-mode gas as having 
Tmax < 10*'^ K. Using this definition, practically all H I ab- 
sorption is due to cold-mode gas. This is not so surprising for 
the Lyman-Q forest, which is comprised mostly of gas outside 
of haloes and would therefore not be expected to have been 
shock-heated to high temperatures. LLSs are dominated by 
low-mass haloes with Mhaio < lO" Mq and Tvir < lO'^ * K. 
Even if this gas were heated to the virial temperature, we 



would still classify it as cold-mode accretion. However, most 
neutral gas in haloes more massive than 10^^ Mq cannot 
have been heated to the virial temperature, or we would 
have found that the fraction of the Hi with Tmax > 10^'^ K 
exceeds the fraction of the Hi contributed by haloes with 
mass > 10^^ M0, which increases from about 20 to more 
than 95 per cent from A''hi = 10^* to 10^^ cm~^. However, 
we find that the actual fraction of Hi contributed by gas 
with Tmax > 10^'^ K is only 5-30 per cent for this column 
density range. Thus, it is clear that even for high column 
density absorbers, the H I content is dominated by gas that 
has never gone through an accretion shock near the virial 
radius. 

We can also think of cold-mode accretion as gas that 
has already accreted onto haloes, but not yet onto galaxies, 
without having reached temperatures above 10^'^ K, flow- 
ing in with relatively high velocities, and that will accrete 
onto a galaxy in the near future. About 30 per cent of all 
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cm systems satisfy all these criteria. 
Therefore, even this set of highly restrictive criteria leads to 
the conclusion that cold-mode accretion has already been 
observed in the form of high column density H I absorption 
systems. 

However, if we think of cold-mode accretion as accret- 
ing onto galaxies for the first time, in addition to having low 
Tmax, being inside haloes, and rapidly inflowing, the fraction 
of LLSs satisfying this requirement is about a factor of two 
lower, but still around 15 per cent. Cold streams also contain 
galaxies, so the gas that is being recycled could still be ac- 
creting along these streams. If we include previously ejected 
gas in our definition of cold-mode accretion, the metallicity 
of c old-mode gas wo u ld not necessarily be low. 

iFumagalli et al.l |2oi3) already pointed out that cold 
streams have likely been detected as LLSs. They used 
high-resolution simulations of seven galax;ies with 10^"'^ < 
Afhaio < 10^^'^ Mq at z « 3 by resimulating a cosmologi- 
cal simulation. They included ionizin g radiation fro m local 
stellar sources, which we did not. iNagamine et al.l ((20ioj) 
have shown that local stellar sources change the column 
density distribut i on fu nction by less than 0.1 dex, although 
IFumagalli et al.' ("2011' ) find a so mewhat larger effect (0- 
0.5 dex). Fuma galli et al.l l|201ll ) investigated Hi absorp- 
tion from "central galaxies", which they defined as all gas 
within 0.27lvir, and from "streams" , defined as all gas at radii 
0.2i?vir < R < 2RviT, and found that "cold streams" mostly 
have 10" < A^'hi < lO" cm . However, their definition of 
"cold streams" was all gas at 0.2_Rvir < R < 2_Rvir, which 
excludes neither gas that has been shock-heated to the virial 
temperature, i.e. hot acc retion, nor gas that is o utfiowing. 

The halo sample of IFumagalli et al.l (|201lh underpro- 
duces A''hi < 10^^ cm~^ systems, which they argue is be- 
cause many of these absorbers live outside of haloes. This 
is in agreement with our finding that the contribution of 
halo gas is strongly declining with decreasing A^hi, although 
we find that up to 30 per cent still stems from low-mass 
haloes not covered by their simulations. We find that for 
10^^ < iVni < 10^° cm-^ haloes with 10^ < Mhaio < 
10^" M0 account for about as much Hi absorption as all 
higher-mass haloes combined. The fact that 10^* < A'hi < 
10^" cm~^ systems are also underproduced in their sample 
is therefore consistent with our finding that very low-mass 
haloe s dominate the H i colu mn density distribution func- 
tion. iFumagaUi eFai] (|201ll ) conclude that haloes in the 
mass range 10^" < Mhaio < 10^^ M© already account for 
20-30 per cent of the LLSs and DLAs. We have found this 
fraction to be even larger: 30-40 per cent for LLSs and 40- 
70 per cent for DLAs. This difference could be due to their 
inclusion of ionizing radiation from local sources. 

Our results lead to the somewhat stronger conclusion 
that for 10^* < A''hi < 10^^ cm"^ more than 80 per cent 
of the H I absorption is caused by cold- mode halo gas. Only 
about 5 per cent is caused by hot-mode gas. The remainder 
may arise from the high-density part of the IGM or from 
haloes with Mhaio < 10^ Mq. At lower column densities 
(A''hi < 10^* cm~^), the contribution of halo gas declines, 
though more than 95 per cent of the absorbing gas is still 
cold mode. For higher column densities (A''hi > 10^^ cm~^), 
the absorption is dominated by the ISM of galaxies in haloes 
with Mhaio > 10^^ Mq and, depending on the column den- 
sity, the mean cold-mode contribution is between 70 and 



90 per cent. At all column densities, inflowing gas accounts 
for most of the H I column density, but there are significant 
contributions from static and outfiowing gas. 

We conclude that cold streams are real. They have al- 
ready been observed in the form of high column density ab- 
sorbers, mainly in systems with 10^^ < A''hi < 10^^ cm~^. 
Cold fiows are critical for the success of our simulation in 
reproducing the observed z = 3 column density distribution 
of damped Lyman-a and particularly that of Lyman limit 
systems. 
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